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ABSTRACT 

Time-series data taken with the Hubble Space Telescope of a field six core 
radii (~ 5') from the center of the globular cluster Messier 4, and covering a period 
of about 10 weeks in early 2001, have been analyzed in search of variable objects. 
Various criteria were employed to select candidate variable stars. Period searches 
were performed on the selected candidates using phase dispersion minimization 
(PDM). The reliability of the PDM search results was tested using synthetic light 
curves of eclipsing binary stars and sinusoidal light curves of different periods. 
Results from this analysis showed that there are probably no eclipsing binary stars 
or periodic variables in our field with periods on the order of a few hours to a few 
days, down to limiting magnitudes of ~ 25 and / ~ 24, which is consistent with 
the absence of contact binaries such as W Ursae Majoris systems. However, one 
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candidate variable star does show an increase in brightness of ~ 0.1 magnitudes 
in both bandpasses, which seems to last for a few days. Possible explanations 
concerning the nature of this object include a binary system with a white dwarf 
primary and a low-mass main sequence secondary, or a BY Draconis variable 
star. We are able to set an upper limit to the observed fraction of photometric 
variability in this dataset of 0.05%. 

Subject headings: binaries: eclipsing — stars: variables — globular clusters: indi- 
vidual (M4) 

1. Introduction 

It is well known that binary stars are important for the determination of stellar masses. 
In addition, knowledge of the population of these objects in globular clusters (GCs) is funda- 
mental to understanding the evolution and dynamical history of GCs. Even a small fraction 
of binary stars in a cluster can contribute greatly in preventing or delaying core collapse of 
the cluster. A full treatment of this topic can be found in Hut et al. (1992). 

This paper presents the results of a search for variability in faint objects in a time-series 
dataset taken with the Hubble Space Telescope (HST) in a field of the globular cluster 
M4 (NGC 6121; a = 16'^23"^54!61, 5 = -26°32'23'.'93 (J2000)). Even though the the data 
were not primarily collected to search for variable stars, the potential for the discovery of 
such objects did exist, since it included 246 separate observations of the same field over 
about two and a half months. These observations represent the deepest ever imaging of a 
globular cluster, reaching the faintest part of the white dwarf sequence ever observed. In 
the context of this project, it allowed for the study of variability within magnitude ranges 
that have seldom been explored. Another recent study was performed with HST by Albrow 
et al. (2001) in which nearly continuous observations for 8.3 days were taken in the core of 
the globular cluster 47 Tucanae for the primary purpose of searching for planetary transits. 
Although no such transits were found, the study did identify 114 variables, most of which 
are believed to be binary systems. 

In the color-magnitude diagrams of globular clusters, observations of a binary star se- 
quence 0.75 magnitudes brighter than the main sequence are rare. Even in cases where one 
is found, many stars in the sequence are actually optical doubles, due to the effects of crowd- 
ing in globular cluster images. However, binary stars have been found in globular clusters 
(see, e.g., Hodder et al. (1992), Yan and Mateo (1994), Rubenstein and Bailyn (1996), and 
Kaluzny et al. (1997)). In fact, W UMa type eclipsing binaries are among the most common 
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type of short-period variable stars found in GCs (Rucinski 2000). 

In addition to binary systems, the data were probed for other types of variable objects, 
such as pulsating stars, planetary transits, as well as nonperiodic objects such as flare stars, 
and supernovae in background galaxies. The study of these objects is vital to research in 
diverse areas such as stellar populations, cosmology, and the evolution of stars, clusters, and 
galaxies. 

The following section deals with the observations of M4 and the reduction of the data. 
Section 3 discusses what one may expect to find in these data based on previous studies of 
M4, as well as other globular clusters observed with HST. Section 4 describes the process of 
selecting candidate variable stars, and section 5 discusses how these candidates were searched 
for periodicity. Section 6 explains how the reliability of results obtained from the period 
searches were tested. Section 7 describes searches for nonperiodic variable stars, planetary 
transits, and supernovae. Finally, section 8 presents and discusses results of the analysis. 

2. Observations and Data Reduction 

The data are taken from Hubble Space Telescope (HST) observations of M4 in the pro- 
gram GO 8679 (cycle 9), using the Wide Field Planetary Camera 2 (WFPC2). Observations 
were taken over 123 orbits, covering a period of 68.2 days: 1 January - 9 April, 2001, in 
the natural HST photometric filters F606W and F814W (Holtzman et al. 1995), which will 
be simply referred to here as V and /, respectively. The dataset consists of 98 x 1300 s 
exposures in the F606W filter {V band) and 148 x 1300 s exposures in the F814W filter (7 
band), all taken in the same field of M4, about 6 core radii (~ 5') from the cluster center. 

2.1. Preprocessing and Combined-Image Photometry 

The images were preprocessed according to recipes given in Stetson (1987) and Stetson 
et al. (1998). Crowding was not an issue on these images, due to their distance from the 
cluster center, as well as the high resolution of the HST images. A master star list, frame-to- 
frame coordinate transformations, and photometry of the combined image of the field were 
completed as discussed in Richer et al. (2002), which resulted in color-magnitude diagrams 
that reach apparent magnitudes of F ~ 30 and I ~ 28. 
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2.2. Single-Image Photometry 

To create time-series data for each star found in the M4 field, it was necessary to perform 
photometry on each individual exposure. To accomplish this, the following steps were taken: 
(1) the known frame-to- frame transformations were used to create a combined median image 
of the field for each chip on WFPC2. This was performed using the M0NTAGE2 software 
written by Peter Stetson (1994), which transforms the coordinates of each individual image 
to match a reference coordinate system, and then produces a combined image. (2) Once 
the combined median image was created, aperture photometry was performed on all objects 
using DAOPHOT. Subsequently, point spread function (PSF) photometry was performed 
using the ALLSTAR program in order to identify stellar centroid positions. This was done 
by employing PSFs built specifically for each chip and filter on WFPC2 by Stetson (2001). 
(3) The output photometry and frame-to-frame transformations were then provided to ALL- 
FRAME, again written by Stetson, in order to perform PSF photometry on each individual 
exposure (Stetson 1994; Turner 1997). Once this single-image photometry was completed, 
the location of any given star and its accompanying photometry was specified for every single 
exposure taken in the M4 field. 

Since HST is above the Earth's atmosphere, it is not subject to seeing effects experienced 
by its ground-based counterparts. However, it is susceptible to slight focus changes due to the 
telescope moving in and out of the Earth's shadow as it orbits, which causes temperature 
fiuctuations (i.e. breathing). As a result, the full width half-maximum (FWHM) of the 
stellar PSFs on the frames can vary slightly between individual images. The result is that 
small offsets in the magnitudes calculated from performing PSF photometry will exist from 
frame to frame. These must be corrected before performing any analysis on the data. To 
accomplish this, the magnitudes that are output by ALLFRAME for every image on a given 
chip were compared to those of a reference image, which was chosen to be the first image 
taken in the dataset^. The magnitude offset between frames was calculated by taking the 
weighted mean of the magnitude difference for each star common to both frames. 

The usual way of rejecting data points consists of fitting the distribution of magnitudes 
of a given star with a Gaussian profile, and throwing out any data points at a magnitude 
greater than 3a or so from the mean of the distribution, and perhaps reiterating this process 
a number of times. When looking for variability, objects searched may include those which 
vary greatly in magnitude, such as eclipsing binary stars with deep minima, or supernovae. 



^Use of the combined median image as a reference was also explored, but there was a negligible difference 
found in the final magnitudes and magnitude errors obtained, compared with those found through use of a 
single exposure as the reference image. 
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By performing this so-called "/ccr clipping" type of data point rejection, data points that may 
be rejected could actually be caused by real effects due to a significantly varying object. In 
order to avoid the loss of potentially useful data, wc based the rejection of a given data point 
on its magnitude error rather than its magnitude. A "bad data point" is considered to be one 
which has a magnitude error that is significantly larger than the rest of the data points for the 
corresponding star's light curve. To illustrate how "significantly large" was defined. Figure 1 
shows plots of root-mean-square magnitude (RMS or a) vs. mean magnitude in I and V for 
each chip. It can been seen that the scatter of data points begins to increase significantly 
at a well defined "elbow" in the distribution, fainter than which the RMS increases rapidly. 
With a few exceptions, the maximum scatter of the light curve of any star found in the 
dataset is approximately 0.1 magnitudes. To be conservative, this is the value chosen as the 
cutoff magnitude error (as returned by DAOPHOT for a single measurement), so that any 
individual data point with an error above this value was rejected, and any data point with 
an error below this value was retained. 



3. Expectations 

In a past search for binary stars in M4, seven W UMa binary stars thought to be cluster 
members were found by Kaluzny et al. (1997) in an 8'.8 x 8'.8 field centered on the cluster 
with apparent magnitudes reaching as faint as V = 17.7. The field observed in the current 
dataset is well within this area. However, unlike the above study, the magnitude range that 
can be probed for variability in the current study ranges from 20 < < 26, so all variables 
found in that study are saturated in the current dataset. However, there is no reason to 
expect that the existence of binaries is exclusive to bright magnitudes. In a recent study of 
the core of the globular cluster 47 Tucanae (47 Tuc), which probed magnitudes as faint as 
V = 25, 11 detached eclipsing binaries and 15 W UMa systems were discovered in a sample 
of 46422 main sequence stars, resulting in an observed binary frequency of 0.056% (Albrow 
et al. 2001). In this study of M4, 2102 stars have been analyzed. If no assumptions arc 
made about the dependence of the distribution of binary systems with cluster radius, one 
would expect that only one or two binary stars should be discovered in the field observed 
based on the observed frequency in 47 Tuc. This already small number is decreased if it is 
assumed that the majority of binaries are expected to reside close to the cluster center due to 
mass segregation, since the M4 field is located at about 6 core radii from the cluster center. 
The picture becomes even less optimistic considering that only 9 of the 15 W UMa echpsing 
binaries found in the 47 Tuc study are below the main sequence turn-off point. However, 
71 BY Draconis variables, which are thought to be members of binary systems, were also 
discovered, making up 0.15% of the sample. This converts to a slightly more optimistic 
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expectation of 3 or 4 of these objects being discovered in the M4 field. 

As a way of gauging the number of supernova events that are expected to be seen 
in the current dataset, a comparison is made to a recent study of the Hubble Deep Field 
(HDF), in which two supernovae were discovered (Gilliland et al. 1999). As with M4, the 
HDF was observed with the WFPC2 camera on HST, and just as deep, with a limiting 
magnitude of / ~ 28 (Flynn et al. 1996; Richer et al. 2002). However, unlike the HDF, the 
M4 images contain zodiacal light, in addition to many stars, due to the cluster's proximity 
to the Galactic plane and ecliptic. This results in much higher background levels compared 
to the HDF, significantly decreasing the number of observable faint galaxies in the field. As 
will be discussed in section 7.2, 170 galaxies were identified in the M4 images. By contrast, 
~ 2100 galaxies were identified in the HDF (WiUiams et al. 1996). Scaling the discovery of 
2 supernovae in the HDF to the current dataset, it is estimated that ~ 0.08 extragalactic 
SNe would be expected in this set of images. 

There also exists the possibihty of discovering variable objects which have not yet been 
discussed, such as cataclysmic variables (CVs), flare stars, transiting planetary-sized objects, 
or perhaps a new species of variable star. CVs are very unpredictable, making it difficult to 
quantify how many are expected to be seen in the M4 fleld. However, in globular clusters, 
which have a high frequency of stellar interactions, it is expected that an abundant population 
of CVs would be found (Di Stefano and Rappaport 1994; Grindlay et al. 2001; Knigge et al. 
2002; Townsley and Bildsten 2002), and indeed, CVs have been found in globular clusters 
(see, e.g., Shara et el. (1996); Cool et al. (1998); Neill et al. (2002)). A relatively small sample 
of stars may limit this population in the current observations, but the possible existence of 
CVs in this dataset is not discounted. 

The main science driver for the above-mentioned study of 47 Tuc was to define the 
frequency of close-in gas giant planets orbiting main sequence stars (Gilliland et al. 2000). 
This study reported a null result, even though simulations that were carried out clearly 
showed excellent sensitivity to planetary transits. However, the observations were centered 
on the core of 47 Tuc, which may be too hostile an environment for planets or protoplanetary 
disks to survive. The M4 images were taken far from the core of the cluster, where the 
number density of stars is much smaller, and thus may provide a more suitable environment 
for planet formation. On the other hand, most extrasolar planets are found to be orbiting 
stars with relatively high metallicity. M4 is a more metal poor globular cluster than 47 Tuc, 
perhaps indicating that a smaller fraction of planets will be found in the current dataset. 
If the null result from 47 Tuc is any indication, it is expected that in the present study of 
M4, which contains only ~ 4.5% of the number of stars in the 47 Tuc dataset, no planetary 
transits will be found. In addition, planetary transits found by Charbonneau et al. (2000), 
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and very recently by Konacki et al. (2003), show eclipse depths of ~ 1% for close-in Jupiter- 
size planets. For a relatively bright, non- variable star in this dataset, the scatter in its light 
curve is typically ~ 0.014 magnitudes. This is most likely too large to confirm a 1% eclipse 
depth should one exist in the light curve. On a more positive note, however, M4 is known 
to contain a millisecond pulsar with a planetary companion (Lyne et al. 1988; Backer 1993; 
Rasio 1994; Sigurdsson 1995; Thorsett et al. 1999; Sigurdsson et al. 2003). The presence of a 
planet in this unusual system may be suggesting that planets are rather common in at least 
this globular cluster. 

A further consideration is that the stars for which planetary transits have been observed 
are stars of solar-like mass and radius. However, the magnitude range of the stars that we 
have examined for variabihty is 19 < < 26. This corresponds to a range in stellar masses 
of 0.6Mq > > O.1M0. These stars have smaller radii than the Sun, and one would 
expect that for planetary occulations of such stars, the eclipse depth that would be seen in 
the star's light curve would be significantly deeper than for the two transiting planets found 
to date. In fact, once the stellar mass reaches O.IMq, close to 100% eclipse depths could be 
expected for a system than is close to edge-on, assuming such a system contains a planet 
with Jupiter's radius. 



4. Candidate Selection 

After reducing the data and rejecting bad data points as described in section 2.2, vari- 
ability searches were conducted on the light curves of the stars in the dataset. Described in 
the subsections that follow are various criteria that were employed in order to narrow down 
the number of candidate variable stars in this dataset. After satisfying the criteria, any star 
that contained in its light curve less than six data points was immediately rejected. This 
may seem like a rather generous cutoff value, but when searching for possible supernova 
events, six points could be enough to determine if it is indeed a viable candidate. Stars were 
also excluded from the candidate hst if: 

• The image of the star is saturated on the CCD; 

• The star is close enough to the edge of the image so that it can partly appear and 
disappear over time due to the observational dithering pattern; 

• The stellar image is affected by light from a nearby saturated star; or 

• The star is partially or exactly aligned with a diffraction spike of a bright or saturated 
star. 
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AU remaining candidates that were not excluded due to the above criteria are presented 
and included in tables accompanying each subsection. In total, 2102 stars were analyzed for 
variability. 

4.1. Statistical Outliers 

A simple identification was first made of outliers in standard deviation (RMS magnitude) 
as a function of magnitude. The light curves of these stars have a larger scatter than is 

expected for stars of similar magnitude, and are thus potential variable star candidates. 
These RMS magnitude outlier stars, not including those rejected due to a lack of data as 
described in section 4, arc shown in Figure 2, and candidates chosen by using this method 
arc listed in Table 1. In this and other tables, coordinates of the objects listed arc based on 
work on M4 by Ibata et al. (1999), and coordinates for stars not included in that study were 
calculated using the CCMAP and CCTRAN tasks on the IRAF analysis package. 

While the above method of choosing variability candidates based simply on the scatter 
of the light curves of stars is a good first pass in the candidate search, it is not terribly robust 
against data corrupted by cosmic rays and warm pixels, which are quite common in HST 
images. A more robust candidate selection criterion that attempts to account for these effects 
is the variability index statistic J, a magnitude-independent measure of the confidence that 
a given star shows real signs of variabihty (Stetson 1996). Another attractive feature of this 
method is that it makes use of data in all filters employed in the dataset when determining 
the variability index. 

Values of J were computed for every star, and those with values of J that were obvious 
outliers in the variability index distribution as a function of magnitude were considered 
variable star candidates. These distributions are shown in Figure 3, where stars selected 
as candidates are plotted as open circles. None of these stars passed the additional criteria 
described in section 4, and so no variable star candidates arose as a result of using this 
method. 



4.2. Main Sequence and White Dwarf Sequence Outliers 

As explained in Section 1, not many binary star sequences are found parallel to the 
main sequence in globular cluster CMDs. Looking at the color-magnitude diagram of M4 
(Figure 4, left; Richer ct al. (2002)), one can see no obvious second sequence in the magnitude 
range observed, other than the expected faint white dwarf (WD) sequence. Nevertheless, 
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there are stars that appear to be outliers to the main sequence, lying above it, or between it 
and the white dwarf sequence. The CMD on the right hand side of Figure 4 is the same as 
that on the left, but stars which obviously do not lie on the MS or WD sequence are plotted 
as open circles. From these stars, period searches were performed on those which were not 
excluded according to the criteria listed in section 4. These remaining candidates are listed 
in Table 2. 

The current datasct also includes white dwarf stars that are coincident with the ZZ Ceti 
instability strip. Figure 5 shows the M4 CMD, with the ZZ Ceti instability strip bordered 
by a black box. Unfortunately, the exposure times of the individual images, which are 
1300 seconds for each single observation, are on the order of the typical ZZ Ceti period. 
Any variability in these stars would be "washed out", since the integration time of each 
observation is too long to capture the star at various stages in its pulsation cycle. It is 
useful however to list the stars which are thought to lie within the ZZ Ceti instability strip, 
as future observations can be made specifically to study these objects. Table 3 lists these 
ZZ Ceti variable star candidates as well as their coordinates, magnitudes and colors, and 
Figure 6 presents finder charts for these three candidate variables. 

5. Periodicity Analysis 

Ever since the landmark paper by Stellingwerf (1978) describing the technique, phase 
dispersion minimization (PDM) has been used extensively and successfully in the discovery 
of countless numbers of variable stars and binary systems. Unlike Fourier transform analysis, 
which is particularly well-suited for sinusoidal-type variations, it is not sensitive to the shape 
of a star's light curve. In addition, its implementation is very straightforward. However, 
determination of the statistical significance of results obtained from use of PDM is not so 
simple, as is the case for Fourier techniques. It has been shown that the PDM probability 
distribution does not follow an F distribution (Schwarzenberg-Czerny 1997), as originally 
thought. Methods used here to estimate the rehabihty of findings made by PDM analysis 
will be discussed in the next section. PDM analysis was performed on all stars in the dataset, 
and particular attention was paid to candidate variable stars. 

6. Reliability Tests 

In an effort to understand the reliability of the PDM analysis performed on candidate 
variable stars, a number of artificial stars were created by sampling different types of variable 
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star light curves at the same observation times as the M4 dataset. These artificial stars 
were then added to the WF2 images, and PDM analysis was performed on the light curves 
extracted from the resulting images. 

Sinusoidal light curves were simulated with a range of periods, amplitudes, and mag- 
nitudes. The range of magnitudes simulated were 17 < / < 28 and 19 < V < 28. The 
bright end of these ranges correspond to the saturation limits of the dataset in each band- 
pass. The faint end corresponds to the faintest magnitude reached in the combined / image. 
Six periods were simulated: 0.085, 0.16, 0.42, 2.1, 5.6, and 15.6 days. The shortest period 
simulated corresponds to a value just above the Nyquist frequency for this dataset, and the 
longest period simulated corresponds to approximately the duration of data collection before 
the occurence of the 44.4-day gap in the observations. Amplitudes of 0.002, 0.02, and 0.2 
magnitudes were simulated. In all, 20181 sinusoidal light curves were simulated for each 
filter. 

In addition, two types of eclipsing binary star light curves were simulated by employing 
the widely-used Wilson- Devinney (W-D) code (Wilson and Devinney 1971). Using this 
program, Algol and W UMa-type systems were simulated, each with a 0.42-day period. 
The simulated Algol system has an orbital inclination (i) of 67°, a relative orbital ellipse 
semimajor axis length (a) of 351.6 i?©, and a primary eclipse depth of «i 0.3 magnitudes. 
The simulated W UMa system has i = 47.5°, a = 1.88 and a primary eclipse depth 
of ^ 0.2 magnitudes. A period of 0.42 days was chosen because it exhibits a case of full 
phase coverage for this dataset. In addition, it is within the range of periods typical of close 
binary stars. In all, 1922 binary star light curves were simulated and added to the images 
in each filter. This light curve data was used in conjuction with the ADDSTAR program 
in DAOPHOT (Stetson 2001) to add artificial stars to the actual M4 images at different 
magnitudes. This ensured that the output light curves would have noise levels typical of real 
stars in the image with similar magnitudes. Photometry was performed on these artificial 
stars in the same way as on the original images (see section 2.2). 

For the artificial sinusoids, a recovered period within 1% of the correct input period was 
considered to be a successful recovery. Results from the simulation of the artificial sinusoidal 
light curves are shown in Figure 7. For the curves in both bandpasses with an amplitude 
of 0.2 magnitudes, PDM analysis recovers the input period almost 90% of the time down 
to a magnitude of 24 in / and 26 in V (at which point the recovery rate drops to ~ 64% 
at its lowest). There is no dependence of the recovery rate on input period apart from the 
0.085 and 15.4-day cases. This gives an idea of the periods to which the PDM analysis is 
sensitive at this amplitude given the time sampling of this dataset. For the case of hght 
curves with 0.02 and 0.002-magnitude amplitudes, the ability of PDM analysis to recover 
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the correct input period is significantly lessened. For an amplitude of 0.02 magnitudes, 
the recovery rate in /-band is very poor regardless of period or magnitude. For y-band, 
the period search is also discouraging, although the 0.42 and 2.1-day period recovery rate 
is slightly higher for brighter stars. In the case of the 0.002-magnitude amplitude curves, 
the PDM analysis recovery rate is even less successful; the input periods are generally not 
recovered, regardless of the period or magnitude of the light curve. This is not surprising, 
since the scatter of the fight curves is on the order of this amplitude or greater (see Figure 
1), even for bright stars. 

The results of the tests for the artificial binary star light curves arc found in Tabic 4. 
These tables show the number of stars added at each magnitude, the number recovered, and 
the percentage of those stars successfully recovered in each filter. A recovery was considered 
to be successful in this case if the period found through PDM analysis is within 1% of either 
the correct input period (0.42 days), half the input period, or | the input period. It was 
found that these are the three major groups within which the periods of the artificial binary 
stars were recovered. This is because the light curve folded at half or three halves the input 
period will cause the primary and secondary eclipses to coincide. This caused the PDM 
analysis to occasionally choose these fractions of the input period. For these artificial binary 
star reliability tests, the PDM analysis is highly successful in recovering the input periods 
of the Algol-type light curves down to limiting magnitudes of / ~ 23 and V ~ 24. For 
light curves fainter than this, less than half are successfully recovered. For the W UMa-type 
light curves, these fimits are 7 ~ 24 (~ 65% recovered) and y ~ 25 (~ 89% recovered). 
Below these limiting magnitudes, almost none of the stars were recovered with light curves 
consisting of more than 6 data points (see section 2.2), let alone with the correct input 
period. 



7. Searches for Non-Periodic Variable Stars 

In addition to periodic variables, there may exist objects in the M4 field that vary in 
magnitude, but not periodically, such as cataclysmic variables, supernovae, and fiare stars. 

There is also the possibility of the presence of objects whose period is too long compared 
with the time frame in which the observations were taken (or too short relative to the time 
sampling or exposure times), so that their periodicity may not be evident, as in the case of 
transits of stars by planets. Clearly, tools such as PDM arc not useful in searches for these 
objects. The variability candidates chosen by some of the methods outlined thus far in this 
section (MS/ WD outliers and variabfity index, to be specific) are vafid for these cases and 
those candidates chosen by these methods were searched for the existence of non-periodic 
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variability in the M4 field. 

7.1. Planetary Transits and Median-Smoothed Light Curves 

As an estimate of the typical period lengths for planetary transits that can be expected 
to be seen in this dataset, we take an example of a Jupiter-sized planet. Prom the reliability 
tests described in section 6, we can expect to find periodic transits that have eclipse depths 
of about 20%. To get a 20% dip in the hght curve of such a system, the radius of the parent 
star must be < O.2i?0. This translates into a mass < 0.23Mo (Montalban ct al. 2000), and 
a magnitude V > 23.2 for this dataset. We set a minimum detectable transit time ~ 5600 s 
based on the minimum separation in time between three consecutive exposures in the light 
curve. Wc estimate the orbital speed implied by such a transit by calculating the length 
of the transit to be the time it takes for the planet to pass from second contact to third 
contact. This implies a maximum orbital speed of 31.6 kms~^. Assuming a circular orbit 
(vorb — 27ra/P), the orbital period is then given by: 

P=H^^(^^±iM^ndays. (1) 

'^orb 

This is longer than the length in time over which the dataset was taken, so only one eclipse 
would be expected to be found. This is also a lower hmit to the period, for the following 
reasons. From section 6, the lower magnitude limit at which we successfully recover periods 
over 90% of the time is 1/ ~ 26. This corresponds to ~ O.IMq for main sequence stars. 
A lower mass star has a smaller radius, and so the orbital velocity of the planet must be 
slower in order to have a transit time that can be detected in this dataset, which results 
in a longer orbital period. The same argument holds for a higher orbital inclination, since 
the separation between second and third contacts is smaller. Thus, periodic signatures of 
planetary transits would not likely be found for such systems in this dataset. Searches for 
potential one-time transits in this dataset are discussed below. 

One method used to search for one-time planetary transits and variable star candidates 
with non-periodic variability involved smoothing the light curve of each star found in the 
dataset. This was done by replacing each data point in the light curve with the median of 
that point and its two nearest neighbouring points. This has the effect of reducing the scatter 
in the light curve, so that real variations in the light curve would be more pronounced. In 
addition, median smoothing was performed using four and six nearest neighboring points, 
to search for longer-term variability in fainter stars having noisier light curve data. 

This method was executed as follows. The faintest data point in each median-smoothed 
light curve was selected, and the ratio of the difference in magnitude between that point and 
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the mean magnitude of the light curve to the RMS scatter of the hght curve was calculated. 
Outliers to the distribution of this ratio represent data points that deviate more than is 
expected from the scatter of the smoothed light curve, which would indicate the possibility 
of a real temporary variation in the star's light curve. However, all the outliers found were 
simply due to saturated stellar images, diffraction spikes, or the star being too close to the 
edge of the image, causing it to disappear and reappear in the frame over time. Thus, no 
candidates for planetary transits, or other type of non-periodic variable star, were found 
using this method. 



7.2. Seeirching for Supernovae 

In addition to the median smoothing of light curves, additional techniques were cm- 
ployed to search for possible supernovae (SNe) occuring in the M4 field. Within the dataset, 
there is a 44.4-day gap in the data. For this section, and further discussions pertaining to 
searching for supernova with this dataset, epoch 1 will refer to the data taken before this 
gap, and epoch 2 to the data taken after the gap. The data taken during epoch 1 and epoch 
2 were combined separately for each bandpass using M0NTAGE2 , and PSF photometry 
was performed on each resulting image using DAOPHOT and ALLSTAR in exactly the same 
way as described in section 2.2. The combined epoch 1 images used 100 images in / and 66 
in V , taken over 14 days. The combined image for epoch 2 used 48 images in / and 32 in 
V , taken over 9 days. Two methods employing these combined images were used to search 
for supernova candidates. The first involved comparing the magnitudes of stars between the 
two epochs. Outliers in the distribution of magnitude difference as a function of magnitude, 
shown in Figure 8, were chosen as possible candidates. However, none of the candidates 
resulting from this search passed the exclusion criteria described in section 4. For the second 
method, background galaxies were identified, and the second epoch image was subtracted 
from that of the first epoch (after performing sky subtraction and flux normalization for 
each combined image) in order to search for residual fiux at the position of these identified 
galaxies. There was no evidence for any residual fiux greater than about la above (or be- 
low) the noise levels for the image, except where the nucleus of the galaxy was too bright, 
in which cases the quality of the subtraction was poor. However, some of these bright nuclei 
show stellar-like hght profiles, perhaps indicating the presence of active galactic nuclei. This 
would make these galaxies interesting objects for possible further study. 
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8. Results and Discussion 

After candidates were selected as described in section 4, their light curves were examined 
more closely and were analyzed using PDM to search for periodicity. Of all the variable star 
candidates listed in the tables in the previous sections, only one of these stars (PC, ID #97) 
showed real signs of variability in its light curve. Even for this remaining candidate star, 
indicators of periodicity were not significant enough to be conclusive. It is therefore possible 
that this remaining candidate is non-periodic, or else there is not enough data in the light 
curve of that star to decisively find a period for it through PDM analysis. This amounts to 
an observed photometric variability fraction of 0.05% for this dataset. 

8.1. Remaining Candidate 

There is the still unresolved issue of the remaining candidate star found in the PC field 
of the observations. There is no clear reason that this star should not be considered as 
variable: it is relatively bright, it is not close to any image defects (such as cosmic rays 
or diffraction spikes from saturated stars) and is far from the edge of the image. Figure 9 
shows a combined J-band image of the PC chip observations, in which the star in question 
is circled. In addition, it has been identified as a cluster member by its proper motion. The 
notable feature in the light curve of this star is the apparent rise in amplitude by ~ 0.1 
magnitudes within the first few days of the onset of data collection, as seen in Figure 10. 
This peak is seen to occur in both bandpasses at the same time, further supporting the claim 
that this object is indeed variable. A clear second peak in brightness is not seen, which could 
mean that the star is not periodic. However, it may be that this object has a period such 
that a second peak occured during a gap in the observations, or else would probably have 
to be modulated in amplitude over time in order to explain the lack of additional peaks. 
Inspection of its light curve seems to indicate a possible period in the range of 5 - 9 days, 
but periodicity analysis performed on this star showed no significant result. 

This star has a magnitude on the combined image ofV = 22.6 mag and color V—I = 0.71 
mag, which places it at the bright end and just red of the white dwarf sequence in the CMD 
of M4, as shown in Figure 11. The location of this star in the color- magnitude diagram 
indicates that it may be a binary star, composed of a low-mass main sequence star and 
a white dwarf. In order to explore this possibility further, a simple decomposition of this 
candidate into two stars was performed. We investigated the magnitudes and colors of two 
separate stars, lying on either the main sequence or white dwarf sequence, that would be 
needed to produce a single observed star of magnitude and color similar to this candidate. 
This gives further insight as to the types of stars that might exist in a binary system found 
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in the position of this candidate on the M4 CMD. We find that a binary system consisting 
of a white dwarf with magnitude and color V ~ 22.7, V — I ^ 0.4 and a faint, late-type 
main sequence companion star, with upper limits to its magnitude and color of V > 26.8 
and V — I ^ 3.4, respectively, result in an object with a position on the CMD that is very 
near to that of the candidate star. This is illustrated in Figure 12. The positions of the two 
speculated binary star constituents are shown as the smaller open circles, and the resulting 
binary is shown as a filled circle. A large open circle is centered on the position of the 
actual candidate star. As one can see, the composite star is very close to the position of the 
candidate. Because of the proximity of the candidate variable to the white dwarf cooling 
sequence, it was expected that such an object would be dominated in luminosity by a white 
dwarf. The position of the hypothetical lower main sequence constituent would indicate 
that this object could be a late M-typc star. However, the magnitude and color found for 
this companion are lower limits, since any companion that is fainter on the main sequence 
contributes little to the hypothetical binary system. It is then possible that this companion 
could be an even fainter, lower-mass main sequence star. 

It may be possible that this object is a cataclysmic variable (CV). It is thought that 
such objects can be identified in the color-magnitude diagrams of globular clusters, between 
the main sequence and the white dwarf sequence (Townsley and Bildsten 2002). In addition, 
those that are found close to the white dwarf sequence will be dominated in magnitude 
by the white dwarf component. The light curve of the object does not seem to show the 
characteristics of an eclipsing binary system. However, the possibility does exist that it is 
a face-on binary system that includes a BY Draconis variable. The period range of 0.4-10 
days of the BY Dra stars found in the 47 Tuc study does not discount this as a possibility. 
In addition, that study reported amplitudes for these stars ranging from V ~ 0.001 - 0.061. 
This is not much smaller than the ~ 0.1 magnitude increase seen in the candidate object's 
light curve. A further possibility is that the observed increase in brightness is due to flaring 
on the star's surface, which would not necessarily be periodic, thus explaining the null result 
of the PDM analysis. However, this is unhkely, since the timescale typical of stellar flares 
rising to maximum brightness and diminishing is on the order of hours, whereas the rise in 
amplitude in the light curve of the candidate star seems to last for days. In any case, more 
data would be needed in order to determine the nature of this object's variability. There is 
another star next to this candidate on the CMD. This star (ID also found on the PC 
chip image, was very close to the edge of the frame. Its partial appearance and disappearance 
from the fleld over time, due to the dither pattern of the observations, was responsible for 
the observed variations in its light curve, discounting it as a possible variable star candidate. 
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8.2. A Lack of Variability 

The lack of eclipsing binary stars in the field can be explained in a number of ways. 
Firstly, it may be that any binaries that are present in this field are face-on relative to the 
line of sight, although this possibility is an unlikely one. Another reason may be the fact that 
the observations were taken in a field that is considerably distant from the core of M4. It 
is expected that mass segregration in the cluster will, over time, preferentially cause binary 
systems to sink toward the cluster center, since binary stars are typically more massive on 
average than single stars. This may not have left many binary systems at the large cluster 
radius that was observed in this study. However, a considerable number of ^ 0.6 Mq white 
dwarf stars are seen in the field (Hansen et al. 2002), as are main sequence stars up to 
~ 0.65 M0 (Richer et al. 2002), so it is not unreasonable to expect, for example, binary star 
systems comprised of 2 x 0.3 stars in that same field. Another possibility is that due 
to the high-density environment of globular clusters, low-mass binaries are more prone to 
disruption before they have the opportunity to harden over time through interaction with 
third bodies. Conversely, the relatively large numbers of blue straggler stars in globular 
clusters (30 have been found in M4 (Rucinski 2000)) may represent the coalescence of many 
low-mass binary systems which have acted as a heat source to the cluster in preventing 
core collapse. This could reduce the number of binaries expected in this magnitude range. 
However, the real explanation for the lack of eclipsing binaries in this sample may simply 
be small-number statistics. As speculated in section 3, based on the low fraction of stars in 
the study of 47 Tuc observed to be in binary systems, the predicted number of these objects 
recoverable in this dataset was low to begin with. 

The small number of background galaxies observable in the M4 images is the most 
probable reason that no supernovae were found in the field observed. Previous studies of the 
Hubble Deep Field, whose images are comparable in depth to the M4 field, yielded the dis- 
covery of two supernovae (GiUiland et al. 1999). Due to the substantially larger background 
flux in the M4 observations compared to that of the HDF (see section 3) , considerably fewer 
galaxies were observable, making it unlikely that any supernovae would be found in these 
observations. 

While this research did not produce any conflrmed discoveries of variability in the magni- 
tude range y ~ 19 — 26 for M4, it is far from having been fruitless. Firstly, there is one inter- 
esting candidate variable, for which it seems that more observations will be needed to identify 
its nature. However, we can place an upper limit to the observed fraction of photometric 
variability for this study at 0.05%, within the magnitude (19 ^ ^ 26 and 17 < / < 24), 
amplitude {AV > 0.1 magnitude), and period (0.16 days ^ P ^ 7-0 days) ranges to which 
we were sensitive. 



-17- 



The identification of the white dwarf stars in the M4 field suspected of being ZZ Ceti 
variables may prove to be of considerable value. The discovery of such stars in a globular 
cluster is of the utmost importance to understanding the chemical structure and evolution 
of population II white dwarfs, which currently suffers from a lack of observational data. 
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Fig. 1. — Plots of RMS magnitude as a function of magnitude for stars found in each chip 
of WFPC2 in the M4 images. Each distribution has a well-defined "elbow" above which the 
scatter of data points for the stars rises sharply. The distribution reaches a maximum in 
each case at an RMS of approximately 0.1 magnitudes. The chip and filter are labelled in 
the top left corner of each plot. 
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Fig. 2. — Shown as in Figure 1 (except for rejected stars - see section 4) are the RMS 
magnitude distributions for each chip and filter. Stars chosen as candidate variables are 
labelled as open circles. 
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Fig. 3. — Variability index J as a function of / magnitude for stars found in each chip of 
WFPC2 in the M4 field. Outhers to this distribution are plotted as open circles. 
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Fig. 4. — Left: Proper-motion cleaned color-magnitude diagram of M4. There are no signs 
of a binary star sequence running parallel to the main sequence of the cluster. Right: Main 
sequence and white dwarf sequence outliers chosen as variability candidates are plotted as 
open circles on the same CMD. 
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Fig. 5. — Color-magnitude diagram of M4. The ZZ Ceti instability strip through the white 
dwarf region is indicated by the box. 
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Fig. 6. — Finder charts for the candidate ZZ Ceti stars in this field. Each chart is 7 x 7 
arcsec and the potential variables are circled. 




Fig. 7. — Results of PDM analysis on an artificial sinusoidal light curve for a range of periods 
and amplitudes. A successful recovery was considered to be within 1% of the input period. 
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Fig. 8. — Distribution of magnitude differences between tfie epocfi 1 and epocli 2 combined 
images. Outliers to tliis distribution were cfiosen as candidate supernovae and are marked 
by open circles. 
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Fig. 9. — The combined /-band PC chip image of the M4 field taken with HST (40 x 40 
arcsec) . The remaining candidate from the variable star search is enclosed by a dashed circle. 
It is a relatively isolated star, in that it is not close to a saturated star or a diffraction spike. 
It is also very far from the edge of the image. 
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Fig. 10. — Light curve of the remaining candidate variable star. The first set of rows shows 
the hght curve before the 44.4-day gap in the observations, and the second set shows the 
hght curve after that gap. It can be seen that the hght curve shows a ~ 0.1 magnitude rise 
in its amphtude a few days into the observations in both the / and V filter data, and is a 
highly desirable characteristic in the selection of this star as a candidate. 
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Fig. 11. — The same CMD of M4 shown in Figure 4, with the remaining candidate star 
circled. It is at the bright end and shghtly to the red of the white dwarf sequence of the 
cluster. The slightly more luminous object just to the red of the candidate was very close 
to the edge of the image, and disappeared from and reappeared in the image over time, and 
was thus excluded from being chosen as a candidate variable. 
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Fig. 12. — Color-magnitude diagram of M4, showing the positions of the two stars (smaller 
open circles) needed to produce a binary system that would be observed with similar mag- 
nitude and color as the candidate star (filled circle). The resulting object is very close in 
position to the candidate, which is located at the center of the large open circle. 
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Table 1. RMS outlier candidates for variability search. 



Chip ID RA^ Dec^ Filters'' V ay V-I ay-i Comments 



PC 


1 


16; 


:23; 


;56. 


.26 


-26: 


:32: 


:27. 


.0 


I, V 


22. 


.349 


0. 


.024 


0. 


.831 


0. 


.042 


isolated; very near to 
edge of image 


PC 


2 


16: 


:23; 


:55, 


.28 


-26: 


:32: 


:35. 


.8 


I, V 


23. 


.383 


0. 


.017 


2. 


.360 


0. 


.020 


isolated; close to edge 
of image 


PC 


3 


16; 


:23; 


;54. 


.59 


-26: 


:32: 


:42. 


.2 


I, V 


24. 


.798 


0. 


.011 


2. 


.862 


0. 


.014 


isolated; close to edge 


pp 
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.zo. 


.oo. 


.OD 


-zo. 


.oz. 




1 

, 1 


T/ 
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.OlO 
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.UU 1 
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.oyu 


u. 




isolated 


PC 


132 


16; 


:23; 


:55, 


.10 


-26: 


:32: 


:04. 


.2 


i,v 


25. 


.969 





.033 


3. 


.295 


0. 


.036 


isolated; next to 






































Sciturated star 


WF2 


21 


16; 


:23; 


:53. 


.41 


-26: 


:31: 


:29. 


.1 


V 


24. 


.041 





.034 


2. 


.468 





.036 


next to saturated stax; 
image 


WF2 


24 


16; 


23; 


53 


.56 


-26: 


:31: 


:29. 


.8 


y 


24. 


.952 





.021 


2, 


.806 





.026 




WF2 


73 


16: 


:23: 


:53, 


.87 


-26: 


;31: 


:24. 


1 


V 


24. 


.040 


0, 


.021 


2. 


.552 


0, 


.027 


next to saturated star 


WF2 


120 


16: 


:23: 


:54, 


.35 


-26: 


:31 


:24. 


,6 


I 


26. 


.096 


0, 


.029 


2. 


.418 


0, 


.036 


next to saturated star 


WF2 


277 


16: 


:23: 


:56, 


.44 


-26: 


:31: 


:31. 


,1 


I, V 


24. 


.097 


0, 


.036 


2. 


.839 


0, 


.046 


next to saturated star 


WF2 


395 


16: 


:23; 


:59, 


.37 


-26: 


:31: 


:57. 


.2 


V 


24. 


.923 


0, 


.033 


2. 


.218 


0, 


.035 


close to edge of image 


WF2 


511 


16; 


;24; 


:00. 


.26 


-26: 


:31: 


:48. 


.3 


V 


25. 


.243 





.014 


3. 


.173 





.018 


near to saturated star; 
close to edge of image 


WF2 


515 


16; 


;23; 


;58. 


.25 


-26: 


:31: 


:10. 


.3 


I 


25. 


.377 





.029 


3. 


.046 





.036 


isolated; appears to be 
two overlapping PSFs 


WF2 


585 


16; 


;24; 


:00, 


.68 


-26: 


:31: 


:37. 


.1 


V 


21. 


.197 





.024 


1. 


.368 





.031 


isolated; close to edge 
of image 


WF2 


661 


16; 


;23; 


;56. 


.63 


-26: 


:32: 


:21. 


.7 


V 


25. 


.530 


0. 


.027 


3. 


.006 





.029 


close to bright star 


WF3 


58 


16; 


;23; 


:59. 


.89 


-26: 


:33: 


:21. 


.0 


V 


20. 


.853 


0. 


.011 


1. 


.865 





.012 


isolated 


WF3 


392 


16; 


;24; 


;02. 


.42 


-26: 


:32: 


:59. 


.9 


i,v 


24. 


.117 


0. 


.036 


0. 


.565 


0. 


.039 


isolated; very near to 
edge of image 


WF3 


554 


16; 


;24; 


:00. 


.53 


-26: 


:31: 


:49. 


.8 


V 


21. 


.152 





.015 


1. 


.203 





.020 


isolated 


WF3 


611 


16; 


;24; 


;01, 


.13 


-26: 


:31: 


:49. 


.2 


V 


22. 


.406 


0. 


.025 


1. 


.349 


0. 


.026 


isolated; close to edge 
of image 


WF3 


613 


16; 


;24; 


:01, 


.34 


-26: 


:31: 


:52. 


.8 


V 


25. 


.322 





.023 


3. 


.122 





.025 


isolated; close to edge 
of image 


WF3 


619 


16; 


;24; 


;01, 


.02 


-26: 


:31: 


:45. 


.7 


I 


19. 


.978 





.053 


1. 


.411 





.059 


isolated; close to edge 
of image 


WF4 


11 


16; 


;23; 


:53, 


.74 


-26: 


:32: 


:53. 


.0 


V 


24. 


.614 





.033 


2. 


.711 





.036 


near another star - PSFs 
may be overlapping 


WF4 


53 


16; 


;23; 


;56, 


.65 


-26: 


:32: 


:32. 


.6 


i,v 


24. 


.965 





.021 


2. 


.781 





.022 


isolated; very near to 






































edge of image 


WF4 


55 


16: 


:23: 


:56, 


.58 


-26: 


;32: 


;33. 


,7 


V 


21. 


.333 


0, 


.005 


1. 


.960 


0, 


.006 


isolated 


WF4 


58 


16: 


:23: 


:52, 


.32 


-26: 


;33: 


;13. 


,8 


V 


21. 


.965 


0, 


.024 


2. 


.070 


0, 


.027 


next to saturated star 


WF4 


63 


16: 


:23: 


:56, 


.58 


-26: 


;32: 


;34. 


,8 


I, V 


22. 


.648 


0, 


.009 


1. 


.880 


0, 


.010 


isolated 


WF4 


190 


16: 


:23: 


:53, 


.79 


-26: 


;33: 


;16. 


.9 


V 


24. 


.564 


0, 


.020 


2. 


.638 


0, 


.023 


next to saturated star 


WF4 


345 


16; 


;23; 


:53, 


.61 


-26: 


:33: 


:42. 


.4 


V 


24. 


.964 





.022 


1. 


.436 





.023 


isolated; very near to 
edge of image 


WF4 


427 


16; 


;23; 


;56, 


.32 


-26: 


:33: 


:32. 


.0 


i,v 


25. 


.599 





.016 


2. 


.735 





.020 


near another star - PSFs 






































may be overlapping 


WF4 


550 


16; 


;23; 


:56, 


.72 


-26: 


:33: 


:48. 


.3 


i,v 


24. 


.235 





.030 


0. 


.580 





.032 


isolated; near to 






































saturated star 


WF4 


551 


16; 


;23; 


;59, 


.03 


-26: 


:33: 


:26. 


.9 


i,v 


24. 


.550 


0. 


.013 


2. 


.742 


0. 


.016 


isolated; near to 
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Table 1 — Continued 



Chip ID RA* Dee Filters'' V ay V-I ay-i Comments 



saturated star 



WF4 


555 


16:23:59.14 


-26:33:26.1 


h V 


24.877 


0.024 


2.794 


0.028 


isolated 


WF4 


563 


16:23:59.19 


-26:33:26.5 


I, V 


22.291 


0.031 


1.431 


0.035 


isolated; near to 
saturated star 


WF4 


579 


16:23:55.60 


-26:34:02.3 


V 


22.006 


0.020 


2.118 


0.021 


isolated 


WF4 


581 


16:23:55.30 


-26:34:05.6 


V 


23.336 


0.022 


2.413 


0.029 


isolated; close to 



of image 



''All coordinates are in J2000. 

''Filter(s) in which star was identified as an outlier. 
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Table 2. M4 main sequence and white dwarf sequence outliers chosen as candidates for 

variabihty search. 



Chip 


ID 




RA 




Dec 




V 


av 


V 


- I 


Cv-I 


Comments 


PC 


1 


16: 


;23:56, 


.26 


-26:32:27. 


.0 


22.349 


0.024 


0. 


.831 


0.042 


isolated; very neaj: to edge of 

image 


PC 


97 


16: 


;23:54, 


.33 


-26:32:20. 


.6 


22.639 


0.009 


0. 


.706 


0.015 


isolated 


WF2 


404 


16: 


;23:56 


.79 


-26:31:07. 


.1 


24.480 


0.027 


2. 


.882 


0.029 


isolated 


WF3 


78 


16: 


:23:58 


.95 


-26:32:58. 


.1 


25.262 


0.023 


3. 


.554 


0.027 


isolated; near to saturated star 


WF3 


532 


16: 


:24:02, 


.37 


-26:32:29. 


.9 


24.014 


0.023 


2. 


.786 


0.027 


isolated; next to diffraction 
spike of saturated star 


WF3 


625 


16: 


:24:02, 


.30 


-26:32:08. 


.4 


27.190 


0.027 


2. 


.821 


0.029 


isolated; next to diffraction 
spike of saturated star; very 
near to edge of image 


WF4 


284 


16: 


;23:54, 


.03 


-26:33:29. 


.6 


23.985 


0.018 


2. 


.044 


0.019 


isolated; near two saturated 
stars 


WF4 


327 


16: 


:23:57 


.94 


-26:32:59, 


.9 


23.409 


0.021 


2. 


,572 


0.025 


isolated 


WF4 


335 


16: 


:23:56, 


.22 


-26:33:16. 


.9 


22.937 


0.016 


1. 


.884 


0.018 


isolated 


WF4 


473 


16: 


;23:54, 


.87 


-26:33:53. 


.3 


24.443 


0.010 


3. 


.007 


0.016 


next to diffraction spike of 
saturated star 


WF4 


519 


16: 


;23:57 


.81 


-26:33:31. 


.4 


22.926 


0.013 


2. 


.153 


0.018 


isolated; near to diffraction 
spike of saturated star 
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Table 3. ZZ Ceti variable star candidates. 



Chip ID RA Dec X Y V ay V - 1 ay-i 



WF2 340 16:23:58.88 -26:32:00.0 70.-567 1.56.524 24.515 0.025 0.649 0.028 
WF4 286 16:23:54.13 -26:33:29.0 388.078 660.529 24.621 0.026 0.700 0.028 
WF4 550 16:23:56.72 -26:33:48.3 749.814 482.751 24.235 0.030 0.580 0.032 
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Table 4. Artificial binary star test results. 



% Recovered % Recovered 



/ 


Algol 


W UMa 


V 


Algol 


W UMa 


17.0 


96.0 


94.6 


19.0 


96.0 


95.9 


18.0 


91.9 


91.9 


20.0 


91.9 


91.9 


19.0 


98.6 


98.7 


21.0 


97.3 


98.7 


20.0 


95.9 


97.3 


22.0 


90.7 


97.2 


21.0 


94.6 


95.9 


23.0 


71.7 


93.3 


22.0 


91.9 


93.2 


24.0 


56.9 


90.5 


23.0 


91.9 


90.6 


25.0 


45.8 


89.2 


24.0 


40.5 


64.8 


26.0 


13.4 


32.5 


25.0 


2.8 


2.7 


27.0 


0.0 


2.7 



^Primary eclipse depth is ^ 0.3 mag for Algol- 
type stars and ^ 0.2 mag for W UMa-type stars. 



